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Abstract 

The formation of multi-soliton configurations in Bose-condensed alkali-metal 
atoms is analyzed by using the nonpolynomial Schordinger equation. A train 
of bright solitons is obtained from an axially homogeneous Bose-Einstein con- 
densate by a sudden change of the scattering length from repulsive to attrac- 
tive. We derive an analytical expression for the number of bright solitons 
generated by using this mechanism. The formula generalizes a previous for- 
mula obtained with the ID Gross-Pitaevskii equation. In the second part 
we consider vortex bright solitons, namely cigar-shaped bright solitons with 
a nonzero angular quantum number k along the axial direction. By using a 
variational approach we determine the shape of vortex bright solitons, show- 
ing that the critical number of atoms for the collapse of the vortex soliton 
increases with a larger k. Finally we calculate monopole and quadrupole 
collective oscillations of these vortex bright solitons. 
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I. INTRODUCTION 



In the last years we have investigated various dynamical properties of pure Bose-Einstein 
condensates (BECs) made of alkali-metal atoms at ultra- low temperatures [1]. The recent 
experimental observation of dark [2] and bright [3,4] solitons in Bose-Einstein condensates 
(BECs) has focused our interest on BEC bright solitons. In particular, by using the non- 
polynomial Schrodinger equation (NPSE), an effective equation derived from the 3D Gross- 
Pitaevskii equation (GPE) [5], we have studied the stability properties of bright solitons, 
determining the critical strength beyond which bright solitons collapse [6]. This collapse is 
a direct consequence of the non-constant size of the soliton in the transverse direction, as 
found by solving the 3D GPE and NPSE. Instead, the ID GPE, for which the transverse 
size of the soliton is constant, admits only stable bright solitons [6]. 

In this paper we analyze the formation of bright solitons in a Bose-Einstein condensate 
of 7 Li atoms induced by a sudden change in the sign of the scattering length from positive 
to negative, as reported in a recent experiment [4]. In a previous paper [7] we have show 
that a number of bright solitons is produced and this can be interpreted in terms of the 
modulational instability of the time-dependent macroscopic wave function of the Bose con- 
densate. From the ID GPE we have derived a simple formula for the number of solitons 
that is in good agreement with the numerical results. Here we generalize such a formula 
taking into account the non-constant size of the BEC in the transverse direction. Then we 
investigate the properties of vortex bright solitons, which are bright solitons with a finite 
quantum number of circulation around the longitudinal axis. 

II. FORMATION OF BEC BRIGHT SOLITONS 

At zero temperature the macroscopic wave function of a BEC made of N dilute 7 Li 
atoms can be modeled by the 3D GPE. For a dilute gas the inter-atomic strength depend 
on the single parameter: the s-wave scattering length a s . Note that only for a s < there 
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exist bright solitons, which are shape-invariant waves traveling without spreading due to the 
interplay between the repulsive kinetic energy and the attractive inter-atomic energy 

Following the experiment of Strecker et al. [4], the trapping potential of the BEC is 
modeled by an anisotropic harmonic potential with transverse frequency u± = 2n x 800 Hz 
and longitudinal frequency u z = 2n x 4 Hz. In addition there is a box optical potential that 
initially confines the condensate in the longitudinal direction. 

By using 3D GPE we first calculate the ground-state wave function of the condensate 
with a positive scattering length a s = 100a B , where a B = 0.53A is the Bohr radius. We 
choose a condensate with longitudinal width L = 284.4 /zm and N = 10 4 atoms. Then the 
ground-state wave function is inserted as initial condition for the time-evolution of the 3D 
GPE with a s = —3a B . We observe that a three-body dissipative term has been added to 
the 3D GPE. This term is important mainly during the collapse because it fixes the critical 
density at which the compression ceases (see also [7]). 

In Figure 1 we plot the probability density in the longitudinal direction p{z) of the 
evolving wave-function. The Bose-Einstein condensate shows the formation of 5 peaks. It is 
remarkable to stress that at the first instants the two external peaks become very high but 
the collapse is avoided by the dissipative term that quikly reduces their density. We have 
verified that without the dissipative term the Bose condensate collapses. The dissipative 
term does not alter the number of peaks but only their density. After the formation, the 
peaks start to separate each other showing a repulsive force between them. Eventually each 
peak evolves with small shape oscillations but without appreciable dispersion. 

The formation of these solitonic peaks can be explained as due to MI of the time- 
dependent wave function of the BEC, driven by imaginary Bogoliubov excitations [7]. Ne- 
glecting the dissipative term, in the ID limit of very long cigar-shaped condensate, the 
Bogoliubov longitudinal elementary excitations e& of the axially homogeneous BEC can be 
obtained from the ID GPE. They are given by 




(1) 
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where n = N/L is the linear axial density and g\d = 2(hu±)a s is the ID interaction strength. 
By suddenly changing the scattering length a s to a negative value, the excitations frequencies 
corresponding to k < k c = (167T | (^irfln) 1 ^ 2 become imaginary and, as a result, small pertur- 
bations grow exponentially in time. It is easy to find that the maximum rate of growth is 
at k = k c /y/2. The wavelength of this mode is A = 2n/k and the ratio L/\ gives an 
estimate of the number N s of bright solitons which are generated: 

JN a s L 

N s = ^ — , (2) 

7raj_ 

with a_i_ = {Ti/imuj^)) 1 / 2 . The predicted number N s of solitons is in very good agreement 
with the numerical results of Figure 1 (see also [7]). 

The Eq. (2) is valid in the ID limit but can be extended to a generic 3D case by using 
NPSE instead of ID GPE. In general, the Bogoliubov excitations can be written as 



efc 



where c is the first sound velocity in the longitudinal direction. The sound velocity c satisfies 
the equation 



W-£- (4) 

where /j, is the chemical potential of the Bose condensate. The chemical potential derived 
from ID GPE is /i — g ld n + fuu±, while using NPSE one finds 

= 9idn = huJ± ( 1 \ 
P y/1 + 2a s n 2 \y/l + 2a s n J y ' 

This chemical potential reduces to the ID GPE one in the ID limit, i.e. for 2a s n « 1. 
Moreover, if a s < the chemical potential /i diverges at 2\a s \n = 1, as a consequence of the 
collapse on the condensate [8]. In conclusion, NPSE tells us that the number iV s of bright 
solitons obtained via modulational instability is given by 



where g = 2\a s \/aj_ is the adimensional inter-atomic strength and p = na± is the adimen- 
sional axial density of the condensate. In Figure 2 we plot iV s as a function of gp for a fixed 
value of p. For small values of gp ID GPE and NPSE give the same results but approaching 
gp = 1 the results are strongly different. In fact, NPSE gives a divergent number of solitons 
at gp = 1. As previously stated, only NPSE takes into account the non-constant width of 
the BEC in the transverse direction. Such a width shrinks to zero for a large attractive 
inter-atomic interaction. 



III. VORTEX BRIGHT SOLITONS 

The 3D GPE of a Bose-Einstein condensate of N atoms confined in the transverse direc- 
tion by a harmonic potential of frequency u± and with an attractive inter-atomic interaction 
(a s < 0) can be derived by the following scaled GPE Lagrangian 



^ t + \v 2 -\^ 2 + y 2 )-^\ 2 



1> (7) 



where tp is the macroscopic wave function of the condensate and 7 = Na s /a±, with a±_ 



(h/(mcu±)) 1 / 2 . Note that in the GPE Lagrangian lengths are in units a±, time in units 



_L > 



action in units h and energy in units hu±. 

Using the cylindrical coordinates r = (x 2 + y 2 ) 1 ^ 2 , 9 = arctg(y/x) and z, we write the 
trial wave function of the vortex bright soliton as 

V(r, z, 9, t) = r k e lke exp {-^p " ^ + + ^W^ 2 } , (8) 

where k is the vortex quantum number and a, i], a and (5 are the time-dependent variational 
parameters. In [9] we have used this trial wave function to study vortex states with positive 
scattering length in a harmonic potential. Here we consider a negative scattering length. 
Moreover, in this case the (harmonic) confinement is only in the transverse direction. By 
inserting the trial wave function in the GPE Lagrangian, after spatial integration, one gets an 
effective Lagrangian depending on the variational parameters. By using the Euler-Lagrange 
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equations one realizes that the time dependence of the phase parameters a and (3 is fully 
determined by that of a and r\ (see also [9]). Thus, one obtains a new effective Lagrangian 
depending on a and r\ only. It is given by 



In the case k = this effective Lagrangian reduces to the one we have successfully used to 
analyze the collective oscillations of the bright soliton in [6]. In fact, the stationary solutions 
of the Euler-Lagrange equations of L e // give the widths of the bright soliton (k = 0) and 
the widths of the vortex bright soliton (k ^ 0). 

In Figure 3 we plot the transverse radial width a and the longitudinal axial width rj 
of the vortex bright soliton as a function of the inter-atomic strength 7. For values of 7 
close to zero the axial width r\ is very large with respect to the radial width a: the soliton 
is a strongly elongated cigar. By decreasing 7 the soliton anisotropy is reduced up to the 
collapse, i.e. the critical point at which the width parameters are no more minima of the 
GPE potential energy. It is interesting to observe that the critical 7 changes with the vortex 
quantum number k: 7 C = —0.78 for k = 0, 7 C = —2.20 for k = 1, and 7 C = —3.59 for k = 2. 
This means that it is possible to avoid the collapse by creating vortex bright solitons with 
a large vortex quantum number k (see also the recent numerical results [10] of Adhikari). 
Another remarkable effect shown in Figure 3 for k 7^ is the crossing of the widths by 
decreasing 7: the vortex soliton becomes slightly disc-shaped. 

In the lower part of Figure 3 we plot the frequencies Q± and S7 2 of the collective oscilla- 
tions of the vortex bright soliton along the axial and transverse directions. In cigar-shaped 
condensates Qi is called monopole frequency and Q2 is called quadrupole frequency [9]. 
These frequencies are obtained from the effective Lagrangian L e ff linearizing around the 
minima of the potential energy (equilibrium solution). The transverse frequency f2 2 remains 
practically constant while the axial frequency fli grows by decreasing 7 but goes to zero at 
the critical point, beyond which the frequencies are complex and the equilibrium solution 
becomes unstable. Note that for a fixed 7, the value of Q\ strongly depends on k. 



L eff = (k + 1)<t 2 + l -r? ~ (k + l)a 2 



a 2 



1 




2 (2k)\ 7 
n2 2k {k\) 2 a 2 V ' 



(9) 
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IV. CONCLUSIONS 



In this paper we have shown that a train of bright solitons can be created in a Bose- 
Einstein condensate by a sudden change of the scattering length from positive to negative. 
The number of bright solitons can be predicted by using an analytical formula we have de- 
rived from the nonpolynomial Schrodinger equation. We have also investigated vortex bright 
solitons, which are bright solitons with a finite angular momentum along the longitudinal 
axis. Our results suggest that vortex bright solitons can be created with a larger number of 
atoms than bright solitons without vorticity. Moreover, our calculations suggest that vortex 
bright solitons can be experimentally detected by analyzing their collective oscillations. 
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FIGURE CAPTIONS 



Figure 1: Axial density profile p(z) of the Bose-Einstein condensate made of 10 4 7 Li 
atoms obtained by solving 3D GPE. For t < the scattering length is a s = lOOa^, while 
for t > we set a s = —3a B with a# the Bohr radius. Length z in units a z = (h/mtOz) 1 ^ 2 , 
density p in units l/a z and time t in units uj~ x . Adapted from [7]. 

Figure 2: Number N s of bright solitons created via modulational instability, p = Na±/ L 
is the initial (adimensional) axial density of the axially homogeneous Bose condensate, g = 
2\a s \/a± is the inter-atomic strength with a± = (h/muj±) 1 ^ 2 and a s < the final scattering 
length. 

Figure 3: Top: transverse width a and axial width 77 of the vortex bright soliton as 
a function of the interaction strength 7 = Na s /a±. Bottom: monopole frequency Q± and 
quadrupole frequency Sl 2 of the vortex bright soliton as a function of the interaction strength 
7. Lengths in units a± and frequencies in units u±. k is the vortex quantum number. 
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